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Abstract 

We derive general results for the mass shift of bound states with angular momentum £ > 1 in 
a periodic cubic box in two and three spatial dimensions. Our results have applications to lattice 
simulations of hadronic molecules, halo nuclei, and Feshbach molecules. The sign of the mass shift 
can be related to the symmetry properties of the state under consideration. We verify our analytical 
results with explicit numerical calculations. Moreover, we comment on the relations connecting the 
effective range parameter, the binding momentum of a given state and the asymptotic normalization 
coefficient of the corresponding wave function. We give explicit expressions for this relation in the 
shallow binding limit. 
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I. INTRODUCTION 



Lattice simulations are used in many areas of quantum physics, ranging from nuclear 
and particle physics to atomic and condensed matter physics [IH3]. In such calculations 
the system is solved numerically using a discrete spacetime lattice over a finite volume. In 
practice this finite volume is usually taken to be a cubic box with periodic boundaries. When 
simulating composite objects such as bound states, the boundaries of the periodic box will 
modify quantum wave functions. This leads to finite volume shifts in the binding energies, 
and detailed knowledge of these finite volume effects is necessary to improve high precision 
lattice calculations. 

In Ref. [1], Liischer derived a formula for the finite volume mass shift of S-wave bound 
states of two particles with reduced mass /i interacting via a potential with finite range R. 
When a bound state with energy —Eb is put in a periodic cubic box of length L, its energy 
in the rest frame is shifted by an amount 

— kL 

Am5 = -3|7r^ + 0(e-^''^), (1) 



where k = \/2hEb is the binding momentum and 7 is the asymptotic wave function normal- 
ization defined hj ipsir) =7 e~'^^ / (-\/47rr) for r > R. For potentials with exponential fall-off, 
V{r) ~ exp(— r/i?) for large r, the formula is modified by exponentially small corrections 
when the binding momentum k, is smaller than 1/R. Bound states moving in a finite peri- 
odic volume also have a topological phase correction to the energy [HI E] . This topological 
phase correction contains information about the number and mass of the constituents of the 
bound states. It must be included when determining scattering phase shifts for composite 
objects in a finite volume. 

In Ref. [7j, we briefly discussed the generalization of Liischer's formula ([T]) for the finite 
volume mass shift to higher partial waves. In this paper, we present the full derivation of 
these results. Moreover, we give the general mass shift for states with angular momenta 
up to £ = 3. In general, the mass shift for a given state depends on its transformation 
properties with respect to the symmetry group of the cubic box. In addition to reducing 
finite volume effects from precision lattice calculations, our finite volume results can also be 
used as a diagnostic tool to probe the angular momentum and radial structure of the bound 
state wave function. A summary of our results is given by Eq. (66) together with Table [l] 



at the end of Sec. IV A We also derive a general mass shift formula for two-dimensional 
systems in a finite with periodic boundaries. 

In our discussion here we only consider stable bound states. However, resonances above 
continuum thresholds can also be described from the volume dependence of avoided level 
crossings [H [9]. Our results are universal and can be applied to a wide range of systems. 
We will discuss a few examples below. 

In particle physics, there is some interest in hadronic molecules with angular momen- 
tum [TUHT^ . In the case of S-waves, the deuteron and some exotic weakly bound states 
such as the H-dibaryon were recently studied in lattice QCD [13j. Similar investigations 
for exotic bound states with angular momentum appear feasible in the future. In atomic 
physics, several experiments have investigated strongly-interacting P-wave Feshbach reso- 
nances in ®Li and "^"K [HHIS], which can be tuned to produce bound P-wave dimers. If 
such systems are simulated in a finite volume, our results can be used to describe the finite 
volume dependence of the dimers. 
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Halo nuclei constitute a special class of weakly-bound nuclei with molecular character. 
They consist of a compact core and one or more nucleons with low separation energy [T7] . 
Among halo nuclei there are some systems with nonzero orbital angular momentum. A 
well-known example of a P-wave halo state is the = l/2~ excited state in ^^Be. The 
electromagnetic properties of the low-lying states in this nucleus can be well described in a 
two-body halo picture of a ^°Be core and a neutron [181 HH] • A similar strategy was applied 
to radiative neutron capture on ^Li [20] . A related class of systems is given by nuclei with 
an a-cluster structure such as ^Be and excited states of ^^C [2TVI23] . The Hoyle state in ^^C 
was recently calculated for the first time in an ab initio approach using chiral effective field 
theory and nuclear lattice calculations [23]. Understanding the finite volume corrections 
is crucial in such a calculation. Due to the a-cluster structure, the volume dependence of 
the Hoyle state is governed by the volume dependence of a three-body system. Finally, we 
note that the asymptotic normalization coefficient of the bound state wave function appears 
in the mass shift formula. Hence, our results can be used to extract this quantity from 
lattice calculations at finite volume. The asymptotic normalization is directly connected to 
zero-energy capture reactions [21], which play an important role in nuclear astrophysics. 

The paper is organized as follows. We start with some prerequisites and a general discus- 
sion of the finite volume mass shift in Section [Tll Liischer's result for S-waves is reviewed in 
Section HI, while our extension to higher partial waves is given in Section IV In particular, 
we discuss the mass shift for the irreducible representations of the cubic group, relate the 
sign of the shift to the leading parity, and derive a trace formula for the multiplet-averaged 
mass shift for states with arbitrary angular momentum i. In Section |V| we verify our results 
numerically for two model systems. The case of two spatial dimensions is treated in Sec- 
tion VI Finally, a brief summary and outlook is given in Section VII In the appendix, we 
discuss the relation of the asymptotic normalization constant to scattering parameters and 
derive explicit shallow binding relations between the effective range, the binding momentum 
and the asymptotic normalization coefficient. 



II. BOUND STATES IN A FINITE VOLUME 

As a starting point, we first review several results from [1]. We closely follow Liischer's 
derivation, but consider a more general system with arbitrary angular momentum and non- 
local interactions. 



A. Definitions and basic identities 

We start with some definitions and basic identities about bound states in rotationally 
symmetric potentials V{r). This will also help to define our notation used throughout the 
discussion. For a bound state with energy E = —Eb and angular- momentum quantum 
numbers {i,m), we separate the wave function ipi^) as a product of the radial wave function 
M^(r) and spherical harmonics Yf^{6,(f)), 

^(r) = Reir)Yr{9, <P) = "^¥^{9, ^) . (2) 



3 



This leads to the radial Schrodinger equation 

UfXr)=0, (3) 



^ ' - U{r) + 



with p = y/2fiE = ^y—2flEB = Ik and U{r) = 2fiV{r). The mass is denoted by /i since 
in later applications it will be the reduced mass for a two-particle system. The normalized 
radial wave function satisfies 

dr \ui{r)\'^ = 1 . (4) 



For vanishing potential, we have two linearly independent solutions of the free wave 
equation. These are the Riccati-Bessel functions je{pr) and Riccati-Neumann functions 
fiilpr), which can be expressed in terms of ordinary Bessel functions. It is often more 
convenient to work with the Riccati-Hankel functions 

hf{z) = hi{z) ± ije{z). (5) 

These have the asymptotic form e^'^ as \z\ — > oo. For future reference, we give the explicit 
expressions for h^{z) for i = 0, 1, 2: 

(6a) 

e'(^-"/') , (6b) 
4)e^(-). (6c) 



Asymptotic form of bound state wave functions 

In order to discuss the asymptotic form of bound state wave functions we define solutions 
X^p(r) of ([3]) that fulfill the condition 

Xtir)^hf{pr). (7) 
The bound-state solution normalized according to Q can then be written as 

Ui{r) = i^-f Xtp{r) , (8) 

where p = in and 

/ roo \ -1/2 

7= / dr|xt(r)n (9) 





(^) 














ht 














ht 
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(^) 




z 



is the asymptotic normalization coefficient.^ The factor i^ adjusts the phase such that ug{r) 
is a real function. If the potential has a finite range R, i.e., U{r) = for r > R, we have 
the exact identity 

ue{r) = 1^7/1/ (i/tr) for r > R. (10) 



^ For more details, see Ref. [25j, Chaps. 11 and 12. 
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The asymptotic normalization 7 is itself an interesting quantity. In low-energy astropliys- 
ical reactions it determines the capture rate of a single proton or neutron at zero relative 
energy [21]. In the shallow binding limit it is also directly related to low-energy scattering 
parameters. For £ = 0, 

2 1 



whereas for angular momenta 



7^ K 

> 1 one has 



as 



K 



+ 



2k 



2£ 



0{k) as K — )■ . 



;iia) 



(lib) 



We note that 7 will in general depend upon the binding momentum k. We will derive 
the relations (11a) and (lib) in Appendix |Aj The S-wave formula is well known and given 
in although without the 0{k) correction estimate. For higher angular momentum the 
situation is qualitatively different due to the divergence of hf{z) ai z = 0. 



B. Infinite volume 



We now consider a system of two spinless particles with reduced mass fi with zero total 
momentum. We will work with relative coordinates r = ri — r2. The interaction is given by 
a potential l^(r,r'), which we in general allow to be non-local in configuration space. It is 
assumed to be Hermitian, rotationally symmetric, and to have a finite range R, i.e., 



V{r,v') =0 if \r\> R or |r'| > R. 



(12) 



We consider the case when the system has a bound state lips) with angular quantum 
numbers {i,m). The Schrodinger equation, 



H\ijB) = -EB\ijB) 



can be written as 



l^Ar^B{r) + j dVr(r,r')^B(rO = -^iJ^B(r) 



2/i 

in configuration space. We note that for a local potential, 

V{y,y') = V{r) 5(3)(r-r ) 



Eq. ( 14 ) reduces to the familiar form 

1 



2/x 



A, + V(r) 



(13) 



(14) 



(15) 



(16) 



Regardless of the locality of the interaction, the wave function ipsi^) has the asymptotic 
form 



^B(r) = iSFr(r/r)^^^^ {r > R) 



(17) 



where k, = \/2fiEB is the binding momentum. 
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C. Finite volume 



We now consider what happens when the two-body system is put into a cubic periodic 
box with length L ^ R. For this problem it is convenient to define a periodic extension of 
the potential 

VL{r, r) =^V{r + nL, r' + nL) . (18) 

neZ3 

We take {ip) to be an exact periodic solution of the finite volume Schrodinger equation, 

HLm = -EB{L)\^) , (19) 

with the finite-volume Hamiltonian = Hq + Vl. It is clear that Eb{L) approaches the 
infinite volume eigenvalue Eb and lip) — )■ \iIjb) as L — t- oo. 

We now derive a formula for the finite volume mass (energy) shift, 

ArriB = Eb{oo) - Eb{L) . (20) 

Let us define a state {ipo) by adding together periodic copies of the infinite volume wave 
function in (14), 

{v\^o) = M^) = J2'^Bir + nL). (21) 

n 

This clearly satisfies the periodicity condition. Acting upon this state with the finite volume 
Hamiltonian, we get 

HlM^) = H^^i^Bi^ + n'L) + J2J2f + ^' + ^L)'pBir' + n'L) 

n' n' n 

= Y^i^HoijB{r + nL) + j dV 1/(r + n L, r' + n L) ^^(r + n L) 

+ 5Z / + ^' + + "'^) f 

n^n' J 

= -Eb{oo) ^ ^B(r + n'L) + J] / + ^' + ^^(^ + ^'^) ' 

n' n' n^n' 

(22) 

The final result can be written as 

HL\^,) = -EB{^)\^o) + \ri) , (23) 

where we have defined |?7) as 

n' n^n' 

With the substitution r' — t- r' — nL for each term in the sum, this can be rewritten as 

n' n^n' 
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Due to the finite range of tfie potential we only get contributions from the domain |r'| < R. 
We note that |r' + (n' — n)L| > R when n 7^ n' and R L. Therefore we can use 
the asymptotic form of the wave function and find that \ri) = 0(e~'^^Y This means that 
\iPq) is an approximate solution of the finite volume Schrodinger equation (19) for large L. 
Motivated by this, we write the exact finite volume solution explicitly as 



a\^o) + \^') with 1^') = ^^(e""^) • 



(26) 



We take {ip) to be unit-normalized per volume L^. The same is true of \iPq) up to corrections 
of order e~'^^. We will choose a such that 



(27) 



Consider now the matrix element {ipl Hl li^o)- Acting with Hl on {ipo), we get 

{^\ Hl IV^o) = -Eb{oo) {^\^o) + {^\r]) = -Eb{oo) (V^olV'o) ■ « + (^1^) (28) 



according to (23) and (26). On the other hand, acting with Hl on yields 



{iPl Hl m = -Eb{L) {^\^o) = -Eb{L) (V^olV^o) ■ « • 
Combining these two results we find 

Eb{oo) - Eb{L) = AniB = , . 

We first consider the numerator in this expression. Obviously, 

(V^lr^) = a {i;o\v) + (i^'lv) = « (^olr^) + ^(e-^''^) . 



(29) 



(30) 



(31) 



We note that the factor of a here will cancel the a in the denominator of (30). We can now 
simplify further starting with 

n" n' iiy^n' 

For each n we can make the substitutions r — t- r — nL and r' — )■ r' — nL. These leave the 
integrals invariant, and we get 

i^olv) = J2J2Y1 I / dVV'B(r + (n"-n)L)V^(r,r')^B(r' + (n'-n)L). (33) 



n" n' n^n 

Setting m = n' — n and m' = n" — n yields 



(^ol^?) = C ■ X X j dV j dV i)*B{r + ml) F(r, r') ^^(r + mL) 

m' m^O 



(34) 



where C counts the number of repeated periodic copies. The fact that C diverges simply 
reflects the fact that we are working with periodic wave functions with normalization mea- 
sured per volume L^, and C will cancel in the final result. For the integral to be non-zero, 
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both r and r' have to be close to due to the finite range of the potentiaL From the as- 
sumption L ^ R it then follows that all terms with m' 7^ are suppressed by at least a 
factor of e"^'^'^, and we have 



{Hv) = C-J2 J / dV7/'^(r)l^(r,r')^ij(r' + mL) + C(e-2-^). 

The possible nonvanishing values of |m| are 1, v^, 73, ... . We therefore arrive at 
{^o\ri) = C- J2 j "^^^ j dV>^(r)\/(r,r')^B(r +mL) + C(e-^'^^) . 

|m|=l 

For the denominator in ( [30| , an analogous procedure yields 

(^o|V^o) = C-5^ J d'rrBir)Mr + niL) = C- [l + 0(e-'^^)] 



(35) 



(36) 



(37) 



with the same constant C as above. Combining (36) and (37), the constant C cancels and 
we get 

AniB =J2 J '^^'^ J d^^>B(r) V^(r,r') ^/'^(r' + nL) + 0{e-^^^) , (38) 

|n|=l 

where we have renamed m back to n. 



Eq. (38) is a general result valid for any angular momentum. The dependence of the 
mass shift upon quantum numbers {i,m) will emerge from the wave function ipB and the 
resulting overlap integrals in (38). In the following, we explore this dependence in detail 
and denote the mass shift as Am^i'"^\ 



III. LUSCHER'S RESULT FOR S- WAVES 



For i = the asymptotic wave function (|17|) is given as 

1 Uo{r) 



i^sir) = ^B(|r|) 




An r 



(39a) 



with 



'"o(^) = 7/io (i^r) = 7e for r > R. (39b) 
Due to the finite range R <^ L of the potential we only have contributions with |r' + nL| > R 



in (38). Hence we can insert the asymptotic form for ^pB{r' + nL) and get 

g-K|r'+nL| 



Am 



(0,0) 
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^ /" dV /" dV>^(|r|) V{r,r' 



0{e-^^^) . (40) 



nLl 



We can furthermore use the Schrodinger equation (14) to eliminate the potential. Doing 
this and then renaming r' — > r, we get 



Am 



(0,0) 



|n|=l 



A, 
2/i 



E 



B 



-K|r+nL| 



7 



(41) 



E / dVV^* (|r - nL|)^ [A. - ^ + 0(e-^-) . 
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In the second line we have shifted the integration variable and used partial integration to let 
the Laplacian act on exp(— Kr)/r. Finally, we use the fact that exp(— fi;r)/(47rr) is a Green's 
function for the operator — k^, 

[A.-.^]^ = -5(3)(r). (42) 
This allows us to perform the integral and arrive at 



^^(0,0) ^_V]nj2 rB{\nL\) + 0(e-v^«^) 

^ In|=i (43) 

In the last step we have inserted the asymptotic form of the wave function for ?/'^(|nL|) = 
■ip%{L), and the sum yields a factor of six. This is just Liischer's result ([T]) as given in the 
introduction. 



IV. EXTENSION TO HIGHER PARTIAL WAVES 



We now discuss the generalization of the mass shift formula to arbitrary angular momen- 
tum. The general form for the asymptotic wave function is 



r'yh'j{iKr) 



(44) 



S-wave case, we find 



Inserting this into (|38|) and performing steps analogous to those presented above for the 

1^7/1^ (i/tr) 



|n|=l 



2/i 



[Ar-K']rB{r-nL)\Y, 



+ C(e-^"^). (45) 



The crucial ingredient is the relation 





-g-Kr- 








r 



where are the solid harmonics defined via W['{x,y^z) = -R™(r) = r^Y^™( 
derivation of (46) follows from Lemma B.l in Ref. [2H], which proves that 



RTi^)f{r) = RTir){~] fir) 



for any smooth function /(r). We obtain (46) by using the relation^, 

Id" 



z dz 



h^\z) = {-lY z^' hf\z) 



(46) 
The 

(47) 
(48) 



^ The relation (48) is just a special case of Eq. (10.1.24) in [57], which also holds for other spherical Bessel 
functions. 
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and noting that e '^^ = /i[j'(iKr) and h^{z) = izhf'\z), where h^'iz) is a spherical Hankel 
function of the first kind. 



We can illustrate (46) using an example. For the case £ = 1 and m = we have 

1 



hi{iK,r) 



1 + 



(49) 



and 



oc cos 6. A straightforward calculation shows that indeed 



cos 6* 1 + 



Kr 



ld_ 
K dz 



(50) 



with cos^ = z/r. 

Using (46) to rewrite (45), we get 



1_ 
2/i 



|n| = l 



lip 



K 



(51) 

We now integrate by parts and pass the Laplacian through the differential operator 
i?^(— V^/k). Since the operators both consist of partial derivatives, this is not a prob- 
lem when the wave function is smooth. We assume that this is the case, with the possible 
exception of a measure zero region that can be omitted from the integral. 

The partial integrations give a factor (—1)^. We can now proceed in exactly the same 
way as for S-waves. We perform one more integration by parts so that the Laplacian acts 
on exp(— fi;r)/r. This yields a delta function times a factor of —Air, and the final result is 
then 



Am 



(£,m) 
B 



27r7 



|n|=l ^ ^ 



(52) 



r=0 



For '4'B{em) insert the asymptotic form (44) since it is evaluated in the asymptotic 

region. 



A. Results 

For i = 1, we find the same result for all three P-wave states, 

— nL 

Amg'°) = Am^J'^^) = 3|7p ^ + 0(e-^«^) . (53) 

When compared to the S-wave case, the sign of the P-wave mass shift is reversed while the 
magnitude is the same. Qualitatively, this means that S-wave bound states are more deeply 
bound when put in a finite volume while P-wave bound states are less bound. This will be 
discussed in more detail later. 
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We next discuss the results for i = 2. From (52), we find 



wfiere 





= -15|7| 




= +15|7| 




= -15|7| 



FUx) 



Fi{x) 



.F|(^)+0(e-v^-) 



1 27 2 

- + 3X + y X 



63 



63 



2x + 9a;^ + 21x^ + 21x^ , 

11 9 2 21 o 21 4 

- + -a; + -a; + —a; + —x 
4 2 4 4 4 



(54) 
(55) 
(56) 

(57) 
(58) 
(59) 



We note that the size and even the sign of the mass shift depends on the quantum 
number m. To understand this effect, we need to take into account that our cubic finite 
volume breaks the rotational symmetry group down to a cubic subgroup. 



Representations of the cubic group 

The cubic symmetry group (9 is a finite subgroup of 50(3) with 24 elements. There are 
five irreducible representations of O. They are conventionally called Ai, A2, E, Ti and T2, 
and their dimensionalities are 1, 1, 2, 3 and 3, respectively. Irreducible representations 
of the rotation group 5*0(3) are reducible with respect to O for i > 1. For further details 
about the decomposition see, for example, Ref. [9]. 

In our discussion we assume that the infinite volume system has no partial wave mixing, 
such that orbital angular momentum £ is a good quantum number. We also assume that 
there are no accidental degeneracies in the bound state spectrum, so we can use £ as a label 
for the family of cubic representations split apart at finite volume. Parity invariance remains 
unbroken by the cubic volume, and we have P = (—1)^ just as in the infinite volume case. 
For clarity, however, we will indicate parity explicitly with ± superscripts in the following. 

With our assumptions, an S-wave state in infinite volume will map onto an Af state at 
finite volume. Also a P-wave triplet will map onto the three elements of the represen- 
tation at finite volume. For D-waves, however, the five D-wave states are split into a 
triplet and an doublet, 

= T+ © E+ . (60) 
In the following we use the notation \r,i;a), a = 1, . . . , dim(r), for the basis vectors of 



the irreducible cubic representations. We can rewrite the finite volume mass shift in Eq. (38) 
as 

Amg^'^'"^ = {T,i;a\V ^ f{nL)\r,i;a), (61) 

|n|=l 

where T(x) is the translation operator by displacement x. We can also calculate the matrix 
elements of Attib in the {i,m) basis. In this case there will be off-diagonal matrix elements 
connecting (£, m) and (£, m') when m and m' are equivalent modulo 4. 
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According to Ref. |9], the unitary transformation between the two basis sets for the five 
D-wave states is 



and 



|7^2',2;1> = - 




-1) + |2,1)) 


1^2+, 2; 2) = 




-1)-|2,1)) 


|T2+,2;3> = - 


-71 


-2)-|2,2)) 


|i?+,2;l> = 


|2,0) , 




E+,2;2) = 


71 


-2) + |2,2)) 



(62a) 
(62b) 
(62c) 



(63a) 
(63b) 



So, for example, we have 



A (T-a+.S;!) 



''B 



Am 



(2,1,1) 
B 



-kL 



(64) 



where we have defined 



Am 



B 



|n|=l ^ ^ 



B,{i,m2) 



+ C(e-^"^) (65) 



r=0 



as a straightforward generahzation of ^52^ 



As expected from cubic symmetry, the mass shift is the same for all three states, and 
the same for both states. To summarize our results, we write the mass shift for a state 
belonging to irreducible representation F with angular momentum i as 



Am 



(^,r) 

B 



a 



i^) ■ l7l 



-kL 



(66) 



We list the coefficients ct (^-^) for £ = 0, . . . , 3 in Table |l] 



i 


r 


a{x) 





At 


-3 


1 


Tf 


+3 


2 


-'2 


30x + 135x2 + 315x3 + 315x^ 


2 


E+ 


-1/2 (15 + 90x + 405x2 + 945x3 + 945x^) 


3 




315x2 + 2835x3 + 12285x^ + 28350x^ + 28350x6 


3 


T2 


-1/2 (I05x + 945x2 + 5355x3 _^ ig^^Q^^ + 42525x5 + 42525x6) 


3 




1/2 (14 + 105x + 735x2 _^ 34652;3 + ll340x^ + 23625x5 + 23625x6) 



TABLE I: Coefficient a(x) in the expression for the finite volume mass shifts for 
indicates the corresponding representation of the cubic group. 



o,...,3. r 



12 



B. Sign of the mass shift 



The sign of the finite volume mass shift can be understood in terms of the parity of the 
wave function. In infinite volume the tail of each bound state wave function must vanish 
at infinity. In the finite volume, however, the bound state wave functions with even parity 
along a given axis can remain nonzero everywhere. Only the derivative needs to vanish, and 
the kinetic energy is lowered by broadening the wave function profile. On the other hand, 
a wave function with odd parity along a given axis must change sign across the boundary. 
In this case the wave function profile is compressed and the kinetic energy is increased. We 
have illustrated both cases for a one-dimensional square well potential in Fig. [Tj 



odd 



/ \ 



\ / 
\ / 
\ / 



even 



FIG. 1: Wave functions with even (bottom) and odd parity (top) for a one-dimensional square well 
potential in a box with periodic boundary conditions. The dashed lines give the infinite volume 
solutions for comparison. 



In three dimensions, the situation is slightly more complicated, which can be seen from 
the fact that for i = 2 the sign of the mass shift depend on the representation of the cubic 
group even though the parity is just (—1)^ = -|-1 for all states. In order to understand this, 
we consider the basis polynomials for the cubic representations. These basis polynomials 
are obtained by decomposing the cubic basis vectors in terms of solid harmonics which are 
homogeneous polynomials in x, y and z. For i = 0, ... ,4 the basis polynomials are also 
given explicitly in [2H]- 

For a given polynomial P{x, y, z), we define its leading parity as 

IpP = (-1)'^-^% (67) 

where 

c/max = max{deg^ P, degy P, deg^ P} (68) 

is the maximum degree of P with respect to any one of the three variables. It is this 
leading parity that determines the asymptotic behavior of the mass shift as kL — )■ oo. More 
precisely, we have 

" i^) - (-1)'-^' {^lY~"^'^ as ^ oo (69) 

for the a (^) in Eq. 



66 
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It can easily be checked that this relation holds for all results presented in Table |T} For 
= 2, for example, we have the basis polynomials 



P2,T+ ~ a;?/ , yz , zx, 

r> 2 2 2 2 

P2,E+ ~ a; -y , y - z 



(70a) 
(70b) 



and hence cimax = 1 for the T2 representation and cimax = 2 for the representation. 



C. Trace formula 



for a given i. We can rewrite (52) as 



The expressions for the finite volume mass shift become simpler when we sum over all m 

+ C(e-^'^^). (71) 



Am 



(£,m) 
B 



-1 



-1 27r7 



|n|=l ^ ^ 



r=nL 



Inserting the asymptotic form of the wave function, 



B,(£,m)\ 



r=nL 



i^7 hf{iKr) 



(72) 



r=nL 



and using (|46|) a second time yields 



Am^^'"-^ = (-1)^+1 ■ V ( --V, ) Rr ( --V, 

^ |n|=i V / V 



+ C(e-^''^) . 



r=nL 



Now, from the well-known relation 



m=-t 



rio,<i>)Yri 



2£+l 

4:71 



(73) 



(74) 



and i?^(r) = r^Y, 
carries over to 



we get an analogous expression for the solid harmonics, which then 



t (-iv.) wr (-iv.) /(.) 

m=-l \ / \ / 



^.^(A Yfir) 



(75) 



for any sufficiently smooth function /(r). Finally, we have 



(A, 



k''— (r^^O), 



r 



(76) 



which follows from Eq. (|42|). Putting everything together, we arrive at 

12 

1 



y Am^^ = (-1)^+^ ■ ^ ■ ^ y ^7 (A. 



m=-£ 



|n|=l 



+ C(e"^'^^) 



r=nL 



(77) 



(_l)'^+i(2£ + 1) ■ m + Oie-^^^'^) 

jjiL 
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where the sum just yields a factor of six. Dividing by 2^ + 1, we obtain the average mass 
shift for states with angular momentum £, 



Am 



B 



3|7l 



jj,L 



+ C(e-^'"^) . 



(78) 



Apart from the alternating sign, this average shift is independent of £. 

Eq. (78) can be verified explicitly for the the results presented in Sec. IV A| (c/. Table |l]). 
For £ = 2, for example, one has to average over the three-dimensional representation T2 
and the two-dimensional representation 



V. NUMERICAL TESTS 



In order to verify our predictions numerically, we put the Schrodinger equation ( 19 ) on 



a discrete spatial lattice such that the Hamiltonian becomes an ordinary matrix. We then 
calculate the corresponding energy eigenvalues and eigenvectors. 



A. Lattice discretisation 



We use a hat symbol to denote dimensionless lattice units. For example, we have 

L = L/a and Eb = Ee-a, (79) 

where a denotes the lattice spacing. 

The free lattice Hamiltonian is given by 

3 



fx 



(ri)a(n) - -^r ^ [a\n)a{h + ei) + a\n)a{h - ei)) 
^ 1=1,2,3 



(80) 



where a^(n) and a(n) are creation and annihilation operators for a lattice site n and e/ is a 
unit vector in the /-direction. The corresponding lattice dispersion relation is 



E{q) 



cm) 

2fi 



with the lattice function 



g2(q) = 2 5^ (1 - cosg.) = [l + 0(gf 

i=l,2,3 i=l,2,3 



and the lattice momenta 



q = 27rn/L . 



(81) 



(82) 



(83) 



^ The mapping from the angular momentum eigenstates to the cubic group states is a unitary transforma- 
tion. 
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The binding momentum for a bound state with energy —Eb is determined by 

-(xEb = (1 - cos(-i/t)) = (1 - cosh(/t)) . (84) 



The lattice Green's function for the Hamiltonian (|80|) is 

1 



iqn 



Gih,E) =—y — 



(85) 



We impose periodic boundary conditions by defining the distance f to the origin as 



f (n) 



E 

«=1,2.3 



mm i nf, \ L 



(86) 



B. Methods 



We calculate the mass shift using three different methods: 



1. As a direct difference in energies, Eq. (20), where we use a very large volume (Loo) to 
approximate the infinite volume result. 



2. From the overlap formula (38). 



3. Using discretized versions of (43) and (53), which we obtain by replacing exp(— Kr)/r 
with the lattice Green's function. More precisely, we write the asymptotic bound state 



wave function (17) as 



^^(r) = i^7y^"^(r/r) /i+(mr) e'^" ■ 47rG^(r) (r > R) 
and replace the continuum Green's function 



with the lattice version 



G,(n) =G h 



2fi 



Effectively, this amounts to the replacement 

e--^/L-^47rG'^(L,0,0) 

in the mass shift formula. 



(87) 

(88) 
(89) 

(90) 



The lattice Green's function is also used to calculate the asymptotic normalization 7 from 
the lattice data. This procedure has the advantage of avoiding large lattice discretization 



errors. 
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C. Results 



In the following we report physical quantities in units where the reduced mass is set to 

1. 



Gaussian potential 



We first use a Gaussian potential, 

Kjauss(r) = -K) exp (-rV(2i?2)) 



(91) 



with R = 1 and Vq = 6. This potential does not have a finite range in a strict mathematical 
sense, but the range corrections can be entirely neglected in comparison with other errors 
in our numerical calculation. The smoothness of the Gaussian potential helps to minimize 
lattice discretization artifacts. In Fig. |2] we show the S- and P-wave mass shifts obtained 
with the three methods described in Sec. IV Bl above. The results from the three different 
methods described above agree well for both S- and P-waves. In order to compare the 
dependence on the box size L with the predicted behavior, we have plotted log(L • |Ams|) 
against L (we use the absolute value of Atub since the S-wave mass shift is negative). For 
both S- and P-waves, the expected linear dependence is clearly visible. 
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FIG. 2: S-wave and P-wave mass shifts log(L • lAm^l) as functions of the box size L (in lattice 



units) for a Gaussian potential. We show the results obtained from the direct difference Eq. (20) 



(crosses), evaluation of the overlap integral Eq. (38) (squares), and discretized versions of Eqs. (43), 



( |53[ ) (circles). The dashed lines show linear fits to the overlap integral results. 

When we perform a linear fit to the overlap integral data (dashed lines in Fig. |2]) we 
obtain n = 2.198 ± 0.005, I7I = 11.5 ± 0.2 for the S-wave resuhs and k = 1.501 ± 0.004, 
I7I = 7.0 ± 0.1 for the P-wave results. The values for the asymptotic normalization are 
in good agreement with the results I7I ~ 11.5 (S-wave) and I7I ~ 7.2 (P-wave) that are 
obtained directly from the L^o = 40 data. Inserting the corresponding energy eigenvalues 
into the lattice dispersion relation (84), we find k, ~ 2.211 (S-wave) and k ~ 1.501 (P-wave), 
again in quite good agreement with the fit results. The remaining small discrepancies can be 
attributed to the mixing with higher partial waves induced by the lattice discretization and 
the fact that we have not performed a continuum extrapolation to vanishing lattice spacing. 
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2. Simple step potential 



For a simple step potential, 



step I, 



-Vo9{R-r) 



(92) 



which we use with R = 2 and Vq = 3, the numerical calculation becomes more difficult since 
the discontinuous shape introduces considerable lattice artifacts. Yet we discuss it here 
because it has a strict finite range and we find that for a small lattice spacing of a = 0.2 
the results are satisfactory. In Fig. [3} we show a plot analogous to the one presented for 
the Gaussian potential. Again, the results from the different methods agree well and the 
expected linear behavior is clearly visible. Furthermore, the results from the three methods 
agree well with each other already for smaller L (compared to the results for the Gaussian 
potential), as expected from the fact that the step potential does not have a tail. 
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FIG. 3: S-wave and P-wave mass shifts log(L • lAm^l) as functions of the box size L (in lattice 
units) for a simple step potential. The symbols are as in Fig. [2] 



From fitting to the overlap integral data (dashed lines in Fig. Isl) we obtain n = 2.0636 ± 



0.0005, |7| = 29.17 ±0.06 for the S-wave results and k = 1.6192 ± 0.0009, I7 



for the P-wave results. From the = 40 data we find k ~ 2.0666, I7I 
K ~ 1.6242, I7I ~ 12.8 (P-wave). Given that we do not have error estimates for the L^o 
results, the overall agreement is quite good. 



12.48 ±0.05 
29.6 (S-wave) and 
40 



Finally, we also check our result for the D-wave mass splittings using the step potential 



with a = 0.2. In Fig. we show the mass shift for the D-wave states in both the and the 



representation. Due to the polynomial coefficients a (;^) (see Eq. (66) and Table |l]) one 
does not expect a linear dependence on L for log(L ■ lAm^j). Hence, we simply plot Am^ 
as a function of L directly and do not perform a fit. Nevertheless, we see that (except for 
very small L, where obviously the condition L 3> -R is not satisfied) the agreement between 
the three methods to calculate Am^ is very good and hence conclude that our mass shift 
formula indeed gives the right result also for i = 2. 
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FIG. 4: D- 
of the box 



■wave, mass shift AniB for T2 rep. (left panel) and E"^ rep. (right panel) as a function 
size L (in lattice units) for a simple step potential. The symbols are as in Fig. [2j 



VI. TWO-DIMENSIONAL SYSTEMS 

In this section we derive a formula for the finite- volume (or rather finite-area) mass shift 
of bound states in two-dimensional systems. The results can be used, for example, in lattice 
simulations of cold atomic systems, which can be prepared experimentally to be effectively 
two- dimensional 



In two dimensions, the Schrodinger equation is 

- ^Af'^Bir) + J dV\/(r,r') V'B(r') = -EbM^) 



with 



Id 92 1 92 



(93) 



(94) 



in polar coordinates. States are described by a single angular momentum quantum number 
m = 0, ±1, ±2, . . ., and for the wave function we have the separation 



^i3(r) = Um{r)Ym{0) 



with 



-.imO 



'111 



The two linearly independent solutions of the free radial equation 



d2 1 d 
dr^ r dr 



m 



+ P Umir) = 0, 



(95) 
(96) 

(97) 



are just the Bessel and Neumann functions Jm{pr) and Nm{pr). For a bound state, we have 



—K = —2^Eb, and the wave function has the asymptotic form 



Um{r) = J Kmi^r) for r>R, 



(98) 
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where Km is the modified Bessel function of the second kind. It is related to the Hankel 
function of the first kind, 

H^^\z) = Jm{z) + iNm{z), (99) 



via 



Kmix) 



IX I 



(100) 



As in the three-dimensional case, 7 is the asymptotic normalization. Inserting (100) into (98) 
yields a form which is more similar to the three-dimensional expression, Eq. (10). To make 
the analogy to the calculations in Sec. IV as explicit as possible, we will use the Hankel 
function in the following intermediate steps and only express the final results in terms of 
the modified Bessel function. 

Nearly all of the three-dimensional calculation carries over and we just replace all expo- 
nential terms by Hankel functions. The overlap integral for the mass shift is now 

AmS")= $^y"dVy" dV>^,^(r)F(r,r')^B,m(r' + nL) + o(iif«(v^KL)) . (101) 

|n|=l 



From the asymptotic form of the Hankel function 

~2 



z) ~ 




e'l"" 2 4 j as \z\ ^ 00 . 



TTZ 



(102) 



it is clear that in principle we still have an exponential behavior. In deriving Eq. (101) we 
have used this to write 



O {HO^^i^L)') ~ O (i/W(2i.L)) « O {H0^\V2i.L)) . 



(103) 



In the following, we will simply write the correction terms as 0(e v^'^-^j^ as in the three- 
dimensional case. 



The two-dimensional analogue of the relation (|46|) is 

where Rm{r,6) = r"^Ym{6). This follows from 

^^n^(V2°)/(r) = Rm{r 

and 



2D 



(104) 



r dr 



fir) 



--j iJ«(z) = (-l)-z-™if«(z) 



(105) 



(106) 



The proof for Eq. (105) can be carried out in the same manner as the three-dimensional 
proof in Lemma B.l in [28], where one uses the expansion of e^^'^ (2D vectors) in terms of 
Bessel functions. As the final ingredient, we have 



(107) 
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Using all this in steps completely analogous to those in three dimensions, we get 



Am 



^ |n|=l 



1, 



72D 



(108) 



r=0 



For 171 = (two-dimensional S-waves), this directly yields 



Amg) = -2^Ko(/€L) + C(e-^"^) . 
/i 



(109) 



In fact, Eq. (108) can be simplified further. Inserting the asymptotic form for the wave 
function for ^^^,771 ^^nd using (104) a second time gives 



Am 



(m) 



-1) 



m+l 



|n|=l 



r2D 



( — V 

K 



2D 



i-i/P(i/€r) 



r=nL 



+ C(e-^'^^). (110) 



From (96) and Rm{^) = f"^Ym{6) it is clear that 



27r 



which then yields 



72D 



K 



72D 



-2m 



27r 



111) 



(112) 



for any sufficiently smooth /(r). This is essentially the same relation that we used to derive 
the trace formula in the three-dimensional case, only that here we do not have to sum over 
different m. Together with the two-dimensional analogue of (76), 

{Af)"^ H^^\iKr) = /«2™if«(i«:r) (r ^ 0) , 



(113) 



we then get 



Am 



(m) 



m+l 



\m+l 



^ V — (A^D)' 

2„ ^2m ^ r ) 

^ |n|=l 



2|7l 



i—Hl^\iKr) 



+ 0(e-^'^^) 



r=nL 



(114) 



As we shall see in the following, this is the final result for m = and any odd m, whereas 
for even m 7^ things become slightly more complicated. 

In general, we have to take into account that the finite volume breaks the original planar 
rotational symmetry of the system down to the symmetry group of a square. We find that 
states with the same absolute value of m may mix to form good eigenstates in the finite 
volume. More precisely, we have the symmetric and antisymmetric combinations 



m, ±) = —= ( |m) ± |— m) ) 
V2 
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(115) 



for m 7^ 0. When we calculate the mass shift for these states (in the same way as described 
in Sec. IV A), we get mixing terms of the form 



Am 



(m, mixed) 
B 



-1 



7r7 



|n|=l 



72D 



B,—m \ 



+ 0(e-^"^) . (116) 



r=nL 



Since the condition for the mixing of states is 

2m = mod 4 . 



;ii7) 



they do not play a role for odd m (in fact, they vanish in this case). For even m, however, 
we have to take them into account and find 



as our final result. As an illustration, we give the explicit results for \m\ = 2: 



ft 



Am 



(2,-) 
B 
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(118) 

(119a) 
(119b) 



VII. SUMMARY AND OUTLOOK 

In this work we have derived explicit formulae for the mass shift of P- and higher-wave 
bound states in a finite volume and discussed their decomposition into states transforming 
according to the representations of the cubic group. We have compared our numerical results 
for i < 2 with numerical calculations of the finite-volume dependence for lattice Gaussian 
and step potentials and found good agreement with the predictions. For i > 2, the mass 
shift of a given state {£,m) depends on the angular momentum projection m due to the 
breaking of rotational symmetry. When we average over all m in a multiplet, however, 
the absolute value of the mass shift is even independent of i. The mass shift for states 
in representations of the cubic group, however, is the same for all states. The sign of the 
mass shift can be understood from the leading parity of the representations. Finally, we 
have derived corresponding expressions for the finite volume mass shift in two-dimensional 
systems. 

Lattice calculations also provide a method to extract asymptotic normalization coeffi- 
cients, which are of interest in low-energy astrophysical capture reactions. We have pointed 
out how the asymptotic normalization and binding momentum of a shallow bound state can 
be used to extract the effective range from a simulation. 

Our work provides a general framework for future lattice studies of molecular states 
with angular momentum in systems with short-range interactions. Applications to nuclear 
halo systems and molecular states in atomic and hadronic physics appear promising. An 
important next step would be to include Coulomb effects into the framework in order to 
investigate the volume dependence of bound states of charged particles, which are much 
easier to treat experimentally. In particular, this extension of the formalism is important if 
one wants to describe proton-halo nuclei. Another interesting direction would be to analyze 
the volume dependence of resonances along the lines of Refs. [U |9]. 
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Appendix A: Shallow binding limit 

In this section, we derive the relation between the asymptotic normalization of the bound 
state wave function and the effective range of the corresponding two-particle scattering 
process that we have quoted in Sec. 



It was already pointed out by Liischer in [Ij that the asymptotic normalization 7 of 
the bound state wave function is related to scattering parameters. More precisely, the 
analytically-continued elastic scattering amplitude in forward direction has a pole at the 
bound state energy, and the residue of this pole is proportional to |7p. 

In the limit of shallow bound states, i.e., k — )■ 0, it is possible to make a more direct 
connection to the effective range in the corresponding scattering channel. The crucial in- 
gredients for this are given in a general investigation of Wigner causality bounds [29]. We 
express the elastic scattering amplitude /^(p) in terms of the scattering phase shift Si^p), 

2i 

^'^'^ " p-+Mcoa.(p) - i] ' ^^'^ 

for which we have the well-known effective range expansion 

p''+' cot 6,{p) = -1 + ir,p2 + ■ ■ • , (A2) 

where an and are the scattering and effective range parameters, respectively. A scattering 
state is described by a wave function uf^ir^ that is a solution of the radial Schrodinger 
equation ([3]) for positive real center-of-mass momentum p > 0. As it is done in [29], we 
choose the normalization such that outside the range of the potential (r > i?) we have 

u>f\r) = / [cot 6e{p) Ji{pr) + n^(pr)] . (A3) 

For a bound state, we have p = in and cot 6e{p) = i. Hence, the wave function is 

u['^\r) = {ikY [ije{iKr) + ^^(iKr)] = i^K^h^{iKr) (A4) 

in this case. We immediately see that 



uf\r) = -ue{r) , (A5) 

7 

where M^(r) is the radial bound state wave function from the previous sections (normalized 
to 1). Note that the phase convention is chosen such that the wave functions are real. 
From [21] we have the relation 

r -1 2 

re = be{R) -2 Mm / dr u'f\r) (A6) 
K^o L J 
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with 



2 



TT 



-2£+l 



4 1 /i? 



+ 



2tx 



1 fR 



2£+3 



(A7) 



2r(£-i)r(£ + i) /i? 



TT 



0{a-') . 



For £ = 0, we get 

hQ = 2R + 0{a-^) . (A8) 

Note that a^^ — > as k — )■ 0, and since we only consider finite-range potentials, we can in 
fact write 0{aj^) = 0{k). 



Using the normalization of the wave function we can rewrite ( A6 ) as 

.2e 



be{R) - 2 lim 



dr 



R 



u\ '{r) 



(A9) 



The remaining integral can be expressed as 



dr 



R 



h"' I dr 



(AlO) 



R 



For £ = 0, we get 



dr 



R 



u\ >[r) 



-2kR 



2k 2k 



-R + 0{k) as K^O, 



Together with (A8) this yields 



ro + ^ - - = 0{k) . 

7^ K 



Up to the given order this is equivalent to the relation 

7 



1 - Kro 

from [2n]- For i > 1, the integral is 

2 r(£-i)r(£+i) fR^ 

71 



(All) 



(A12) 



(A13) 



poo 

/ dr 


(") / \ 


IR 





+ 0{k) as K -4 0, 



(A14) 



We see that the leading term exactly cancels the one in (A7) such that we arrive at 

2t^' 



re + 



0{k) 



(A15) 



In all expressions above, the asymptotic normalization of course depends on the binding 
momentum, i.e., •j = only in general this dependence can not be calculated analytically. 
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Using the results that we have derived in this paper, it is possible to determine both k and 
7 from the volume dependence of a given bound state. If this state is sufficiently shallow, 
the formulae derived in this chapter can then be used to get an estimate for the effective 
range of the interaction. 
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